import matplotlib.pyplot as plt 
rocket_M = 101
load_m = 100
div = 1000 
u = 3000
M = rocket_M  # 火箭剩余质量

while(rocket_M < 3000):
    v = 0
    M = rocket_M
    m = (rocket_M - load_m) / div
    print(f"有效载荷:{load_m}\t总质量:{rocket_M:>5}\t比例为:1:{rocket_M/load_m:>5.1f}", end="")
    for i in range(1, div + 1):
        v = v + m / (M - m) * u  # 计算火箭速度
        M = M - m  # 计算剩余质量
    rocket_M = rocket_M + 10
    print(f"\t末速度为:{v:.0f}")
while 1:
    pass